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DEFINITIONS 


d  tank  diameter 

ex,ey  empirical  parameters  to  skew  the  boundary  layer  on  a  rotating  sphere 

g  acceleration  of  gravity 

kx,k.  empirical  parameters  to  model  "free  stream"  velocity  at  wall 

/s  slosh  pendulum  length 

i  time 

x,>',r  tank  fixed  coordinate  system 

C(CV,C,)  dashpot  representing  boundary  layer  viscous  shear  about  *,>■  axes,  respectively 
Cs  dashpot  representing  viscous  damping  of  slosh  modes 

D(Dx.Dy )  dashpot  representing  viscous  damping  of  inertial  wave  modes  along  x,y  axes, 
respectively 

£/,£,  energy  dissipation  rate,  per  nutation  cycle,  for  inertial  and  sloshing  oscillations 

F(FI,Fy,F.)  oscillating  liquid  force  exerted  on  tank  and  its  force  components  along  x,y,z 
axes,  respectively 

lx , Iy , l,  moments  of  inertia  of  rotor 

oscillating  liquid  moment  exerted  on  tank  and  its  moment  components  about 
x,y\z  axes,  respectively 

A/0,A/j  nonslosh  and  slosh  liquid  masses 

Mr  total  liquid  mass  =  M0  +  Ms 

TxJf  moments  exerted  on  tank  wall  by  viscous  boundary  layer 

X0.Z o  location  of  tank  center 

%  Full  tank  filling  percentage 

;5,7. 


moment  of  inertia  ratio  =  spin  axis/transverse  axis  moments  of  inertia 
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a, \h  angular  displacement  of  pendulum  with  complex  amplitudes  a0  and  *4‘o 

ay  tank  angular  acceleration  about  >  axis 

>o  equilibrium  angular  position  of  pendulum 

7)x,ny  angular  displacements  of  rotor 

0  cone  or  nutation  half-angle 

k  nutation  frequency  or  rate 

characteristic  frequencies.  See  Eq.  (A4) 
v  kinematic  viscosity 

p  liquid  density 

u>x,wy  tank  rotation  rates  about  x,y  axes,  respectively 

H  steady  spin  rate  about  Z  axis 
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ENERGY  DISSIPATION  OF  LIQUIDS  IN  NUTATING  SPHERICAL  TANKS 
MEASURED  BY  A  FORCED  MOTION-SPIN  TABLE 


INTRODUCTION 

A  successful  spacecraft  design  —  spin  stabilized  about  the  axis  of  minimum  moment  of  inertia 
and  having  a  large  liquid  propellant  mass  fraction  —  depends  on  long  term  vehicle  control.  Liquid 
motion  in  a  nutating  spacecraft  results  in  kinetic  energy  dissipation  that  increases  coning  motion  in  the 
prolate  spinner.  This  coning  motion,  if  not  controlled,  results  in  a  flat  spin,  or  spin  about  the  axis  of 
maximum  moment  of  inertia.  Coning  motion  must  be  minimized  by  active  nutation  control  to  main¬ 
tain  the  original  attitude.  Consequently,  the  maximum  energy  dissipation  rates  from  the  propellant 
motions  must  be  known  to  size  and  to  compensate  an  active  nutation  control  system  for  these  losses. 

This  report  describes  an  experiment  and  results  of  an  energy  dissipation  study  of  liquid  motion 
using  a  forced  motion  spin  table  Specifically,  this  work  considers  nutating,  smooth-walled  spherical 
propellant  tanks,  without  propellant  management  devices,  which  spin  off  center.  Energy  dissipation 
rates  are  estimated  by  the  resirits  of  tests  having  the  following  objectives: 

•  identify  all  liquid  resonant  frequencies,  and  spin  and  nutation  conditions  under  which  reso¬ 
nances  occur; 

•  measure  oscillating  liquid  forces  and  moments  as  a  function  of  spin  rate,  nutation  rate,  con¬ 
ing  angle,  and  liquid  fill  level, 

•  determine  boundary  layer  shear  at  the  tank  walls;  and 

•  observe  and  quantify  liquid  phenomena  that  might  adversely  affect  prolate  spacecraft  under 
the  control  of  spin  stabilization. 

This  research  also  devises  a  mechanical  model  composed  of  a  pendulum,  a  rotor,  and  viscous 
dashpots  that  can  (a)  simulate  liquid  forces  and  moments,  (b)  predict  liquid  resonances  and  energy  dis¬ 
sipation  rates,  and  (c)  be  used  to  scale-up  test  results  to  flight  conditions  including  gravity  field  com¬ 
pensation.  The  numerical  values  for  the  model  parameters  come  from  the  forced  motion  tests.  An 
equivalent  mechanical  model  of  propellant  motion  is  also  a  practical  tool  in  conducting  simulation  stud¬ 
ies  because  it  can  be  readily  incorporated  in  spacecraft  attitude  control  models. 

The  test  matrix  investigated  a  wide  range  of  variables.  These  include: 

•  spin  rate, 

•  nutation  rate  or  excitation  frequency, 

•  nutation  or  cone  angle, 

•  fill  fraction  or  percent  of  tank  occupied  by  liquid, 

•  liquid  density  and  viscosity,  and 

•  tank  location  with  respect  to  center  of  spin. 

Although  the  vehicle  inertia  ratio  is  not  included  as  a  test  parameter,  the  effective  inertia  ratio  is 
fixed  by  the  spin  rate  and  nutation  rate.  The  tests  actually  investigated  a  wide  range  of  inertia  ratios. 


ZEDD  AND  DODGE 


QUALITATIVE  LIQUID  MODE  DESCRIPTIONS 

A  complete  theory  of  the  kinds  of  motion  that  can  occur  in  a  spinning  spherical  tank  is  not  avail¬ 
able,  but  the  differential  equations  of  motion  suggest  that  two  kinds  of  natural  oscillations  are  possible 
(1,2],  They  are: 

•  sloshing  waves  (free  surface  oscillations),  and 

•  inertial  waves  (internal  liquid  oscillations). 

In  general,  both  kinds  of  waves  produce  oscillating  forces  and  moments  about  the  center  of  an 
arbitrary  shaped  tank.  But  for  a  spherical  tank,  liquid  pressure  can  create  only  a  force;  thus  any 
moment  exerted  about  the  tank  center  can  only  be  due  to  viscous  shear  at  the  wall.  Ordinarily,  viscous 
shear  is  negligible  compared  to  the  effects  of  pressure  [3],  Similarly,  an  oscillating  rotation  of  a  spheri¬ 
cal  tank  about  its  own  center  is  transmitted  to  the  liquid  only  by  viscous  shear  at  the  wall.  For  other 
tank  shapes,  viscous  shear  is  ineffective  compared  to  the  oscillating  wall  motion  normal  to  the  wall  sur¬ 
faces.  One  of  the  items  of  interest  is  to  determine  if  viscous  shear  could  cause  significant  inertial 
waves  in  a  spinning  spherical  tank.  The  following  discussion  of  liquid  motion  is  described  from  the 
point  of  view  of  a  coordinate  system  fixed  to  the  tank  center.  In  this  system,  nutation  causes  oscilla¬ 
tory  translations  along  all  three  axes  and  oscillatory  rotations  about  the  two  axes  that  lie  in  the  plane 
normal  to  the  spin  axis. 

Sloshing  waves  are  characterized  by  oscillations  of  the  free  surface  and  center  of  mass  location  such 
as  to  change  the  potential  energy  of  the  liquid  relative  to  the  effective  gravity  force.  The  effective 
gravity  is  a  vector  combination  of  the  true  gravity  and  centripetal  accelerations.  Sloshing  is  a  dynamic 
interaction  between  effective  gravity  forces  and  inertial  forces.  There  must  be  a  free  surface  and  it 
must  move  up  and  down  through  the  gravity  field.  Mathematically,  sloshing  can  be  analyzed  on  the 
basis  of  an  ideal  liquid  executing  an  irrotational  motion,  that  is,  as  potental  flow.  The  effects  of  viscos¬ 
ity  may  be  considered  later  as  a  boundary  layer  at  the  wall;  the  effect  is  to  provide  some  damping  of  the 
motion,  but  viscosity  does  not  significantly  change  the  slosh  modal  characteristics  or  the  natural  fre¬ 
quencies.  When  one  of  the  lower  frequency  modes  is  driven  at  resonance,  the  damping  of  an  ordinary 
low  viscosity  liquid  is  so  small  that  the  wave  becomes  unstable.  For  an  axisymmetric  tank,  the  instabil¬ 
ity  causes  the  free  surface  wave  to  rotate  around  the  symmetry  axis  but  the  bulk  of  the  liquid  still  has 
an  irrotational  potential  type  of  motion.  A  sloshing  wave  is  primarily  excited  by  unsteady  tank 
translation  — not  rotation  — in  a  spherical  tank. 

When  the  tank  spins,  the  liquid  spins  with  it  after  some  initial  transient  motion  The  liquid 
motion  is  thus  rotational,  and  a  conventional  analysis  would  not  apply.  There  are  some  indications  thai 
a  potential  flow-like-sloshing  can  exist  in  this  rotational  field  and  that  all  the  slosh  resonant  frequencies 
are  greater  than  twice  the  spin  rate.  Assuming  that  sloshing  can  be  created  in  a  spinning  tank,  sloshing 
resonances  can  be  excited  only  if  the  excitation  frequency,  is  more  than  twice  the  spin  rate. 

Inertial  waves  do  not  require  a  free  surface  and  can  occur  in  a  completely  full  tank  (2).  The 
center  of  mass  oscillations  are  small  even  if  the  free  surface  oscillates.  The  resonances  represent  a 
dynamic  interaction  between  Coriolis  forces  and  pressure  forces  in  the  bulk  of  the  liquid  interior.  Iner¬ 
tial  waves  are  circulatory  or  to-and-fro  motions  in  the  liquid  interior,  and  there  may  or  may  not  be  any 
apparent  motions  at  the  free  surface  They  are  excited  by  unsteady  tank  rotations  Inertial  wave 
resonant  frequences,  regardless  of  tank  shape,  are  less  than  twice  the  spin  rate,  fl  (4,5). 

Boundary  layer  shear  is  the  coupling  of  liquid  motion  to  lank  walls  via  the  viscous  liquid  itself. 
This  is  an  energy  transfer  mechanism,  and  the  greater  viscosity  liquids  induce  both  greater  shear  and 
energy  dissipation  rates.  Although  the  boundary  layer  is  a  source  of  damping,  it  may  drive  an  inertial 
wave  to  resonance  in  a  spherical  tank  in  addition  to  driving  bulk  motions. 
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MECHANICAL  MODEL 

To  help  understand  the  test  results  and  data  analysis  described  later,  the  equivalent  mechanical 
model  proposed  for  a  spinning  spherical  tank  is  described  below. 


Figure  1  shows  an  equivalent  mechanical  model  based  on  the  liquid  motion  characteristics  previ¬ 
ously  described.  A  pendulum  is  proposed  to  represent  the  predominant  mode  of  sloshing,  and  a  spin¬ 
ning  rotor  is  proposed  to  represent  the  predominant  mode  of  inertial  oscillations.  (Additional  pendu¬ 
lums  and  rotors  could  be  used  to  represent  additional  modes  of  each  kind  if  necessary.)  For  sloshing, 
the  primary  effect  of  the  steady  spinning  is  expected  to  be  the  body  forces  created  by  the  centripetal 
acceleration.  Since  this  body  force  effect  is  similar  to  that  induced  either  by  a  steady  engine  thrust  or 
by  gravity,  the  pendulum  model  is  similar  in  form  to  that  for  a  nonspinning  spherical  tank  in  a  gravity 
field.  The  oscillations  of  the  pendulum  mass  simulate  the  oscillations  of  the  liquid  center  of  mass.  In  a 
spinning  lank,  however,  the  pendulum  has  two  natural  frequencies,  which  correspond  to  oscillations  in 
the  circumferential  direction  and  in  the  transverse  direction.  For  inertial  wave  oscillations,  the  steady 
spinning  is  crucial  since  such  oscillations  do  not  occur  in  a  nonspinning  tank  (21.  Thus  the  equivalent 
mechanical  element  must  also  spin,  that  is,  it  should  be  in  the  form  of  a  rotor. 


The  natural  frequency  of  the  rotor  is  chosen  by  proper  selection  of  the  inertias  lx>  Iy>  and  I.  to 
duplicate  the  natural  frequency  of  the  inertial  oscillation  mode  of  interest.  A  rotational  viscous  dasitpot 
connects  the  rotor  to  the  tank  to  simulate  the  indirect  excitation  of  the  liquid  caused  by  viscous  shear. 
The  moment  exerted  on  the  rotor  by  the  dashpot  is  not,  however,  a  simple  angular  rale  dependency 
since  it  must  simulate  an  unsteady  Ekman  boundary  layer.  (Details  are  given  in  the  appendix  of  this 
report.) 

FORCED  MOTION-SPIN  TABLE 

The  test  stand  simulates  bodies  whose  motions  can  be  described  by  the  equations  for  rotation 
about  a  point.  Rigid  body  motion  about  a  fixed  point  occurs  when  a  symmetrical  rotor  spins  about  its 
axis  which  itself  rotates  about  a  fixed  point.  Three  angles,  known  as  Euler's  angles,  completely  specify 
the  position  of  the  rotor  The  time  rates  of  change  of  these  angles  specify  the  nutation,  precession,  and 
spin  velocities  of  the  rotor.  The  forced  motion-spin  table  simulates  this  rigid  body  motion,  about  a 
fixed  point,  with  spin  rates  driven  by  the  outer  table,  cone  angle  provided  by  the  gimbal  axle,  and  nuta¬ 
tion  rates  driven  by  the  inner  table  Figure  2  shows  construction  details  of  the  test  stand.  Figure  3 
highlights  tank  motions  and  defines  symbols. 
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OES/GN  PARAMETERS: 

I.  SPIN  KATE;  O  -  407 PPM 

t.  MUTATION  KATE:  40-400  CYCLES  PEP  MIN 
X  NUTATION  ANGLE!  0-3* 


TRUNNION  - 


G/M3AL  AXLE 


SPHERICAL  TANK. 
ACRVC  1C,  12  ID 


STRAIN  0A6FS 


-  NUTATION  DR/VE 
MOTOR 


Fig.  2(b)  —  Construction  details  of  forced  motion-spin  table 
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RADIAL  DISTANCE  FROM  SPIN  CENTER 


The  gimbal  housing  is  mounted  on  top  a  spin  balance  table.  Below  the  table  are  twenty  pairs  of 
slip  rings  used  to  transfer  power  to  and  data  from  the  instrumentation.  During  a  test,  the  gimbal  hous¬ 
ing,  or  outer  table,  spins  about  the  vertical  axis  at  a  constant  rate.  The  housing’s  most  important  func¬ 
tions  are  the  support  of  the  gimbal,  gimbal  axles,  and  overall  fixture  weight. 

The  nutation  drive  wheel  rotates  the  nutation,  or  inner,  table  at  speeds  different  from  the  outer 
table.  Inner  table  spin  direction  is  counter  to  the  outer  table  spin  direction  when  simulating  the  inertia 
ratios  of  prolate  spinners. 

The  lower  end  of  the  beam  support  cone  comes  to  a  point  that  attaches  to  a  linear  actuator  or 
piston-in-cylinder  mechanism.  The  piston  remains  fixed  relative  to  the  cylinder  to  maintain  constant 
nutation,  or  cone  angle.  That  is,  the  piston  pulls  the  lower  end  of  the  beam  support  off  center  and 
holds  this  position.  The  cylinder  base,  in  turn,  spins  with  the  nutation  table  generating  constant  cone 
angle. 

The  center  of  the  gimbal  axle  is  the  fixed  point  center  of  spin.  From  this  point,  radial  arm  length 
and  vertical  height  to  the  geometric  centers  of  each  tank  are  measured.  Tanks  can  be  relocated  along 
the  beam  length  to  vary  the  radial  distance  and  on  top  or  below  the  beam  to  vary  the  height  distance. 

An  analysis  of  the  rigid  body  motion  produced  with  this  forced  motion  spin  table  shows,  for  small 
cone  angles,  the  following  relationship  between  spin  rate,  fi,  and  nutation  rate.  A: 

A  =  -(1  -/,//,)«,  (1) 

hence  to  simulate  a  prolate  spinner  where  Ijl,  <  1,  A  must  be  of  opposite  sense  as  fl  with  a  magni¬ 
tude  given  by  the  above  equation.  For  example,  if  1J I,  —  0.50  and  Cl  =  60  rpm  (clockwise),  then  A 
=  30  rpm  (counter-clockwise).  Figure  4  shows  table  operating  nutation  speeds  that  correspond  to 
specific  inertia  ratios.  Note  the  eight  horizontal  lines  representing  tests  of  constant  spin  rate  and 
counter  nutation  rates. 
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SHOWING  SLOSH  WAVE  RESONANCES 
EXCITED  OUTSIDE  PHYSICALLY 
REALIZABLE  INERTIA  RATIOS 
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Model  Tank  Construction  and  Instrumentation 

Forces  and  moments  imparted  to  the  tank  by  liquid  motions  are  the  fundamental  quantities  to  be 
measured  from  this  experiment.  This  section  reviews  propellant  tank  model  instrumentation  used  to 
accomplish  this  objective.  The  liquid  tank  and  an  equivalent  weighted  control  tank  are  each  supported 
by  a  strain-gage-dynamometer  assembly.  Figure  5  identifies  the  tank  mount  components  and  defines 
the  tank  centered  coordinate  system. 


The  plastic  tank  is  held  between  two  support  rings  by  four  rods  boiled  through  the  corners  of  each 
nng  These  rods  compress  the  rings  around  the  tank.  The  plastic  has  been  slightly  deformed  from 
spherical  shape  reacting  to  the  compression  necessary  to  hold  it  rigid. 

A  clevis  attaches  underneath  each  side  of  the  horizontal  support  ring.  The  top  bracket  is  solid  and 
holds  a  spherical  bearing  previously  pressed-fil  into  its  side.  This  bear.ng  fits  into  the  clevis,  and  a  pin 
is  Iriciion-lit  through  both  clevis  and  bearing  Washers,  used  as  spacers,  prevent  further  sliding.  The 
top  bracket  clamps  to  each  of  two  strain-gage-plates,  which  in  turn  clamp  to  the  bottom  bracket, 
f  inally,  the  four  bottom  brackets  bolt  on  a  circular  plate  that  anchors  to  the  beam. 

The  dynamometers,  each  pair  responding  to  forces  along  one  axis,  are  the  components  between 
the  pin  and  the  arc  dar  plate.  The  reaction  end  of  the  dynamometers  is  pinned  through  the  spherical 
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bearing  to  preveni  thv  introduction  of  end  moments  into  the  dynamometer.  Ideally,  only  unidirectional 
forces  act  on  the  dynamometer,  thus  achieving  pure  tension  or  compression  on  the  sections. 

The  gaging  sections  of  'he  strain-gage  plates  are  only  long  enough  to  accommodate  the  strain 
gages  placed  on  them.  The  remainder  of  the  dynamometer  is  as  rigid  as  possible.  This  combination 
results  in  the  highest  possible  buckling  strength  and  stiffness.  The  system  is  not  responsive  to 
moments  about  the  r-axis;  but  all  other  forces  iFx,Fy,F:)  and  moments  are  resolved.  Five 

strain-gage  bridges  per  tank  sample  the  three  axis  forces  and  two  moment  components  produced  by  the 
liquid  reaction  to  externally  forced  tank  motions.  Linear  accelerometers  mount  on  the  liquid  tank  to 
ptovide  reference  motion  for  calculation  of  liquid  phase  lag.  The  recorded  force,  moment,  and  phase 
lag  magnitudes  vs  the  nutation  frequency  of  excitation  provide  the  basis  for  dissipation  rate  calculations. 
These  recordings  along  with  qualitative  descriptions  of  liquid  motions,  observed  via  video  camera 
mounted  on  the  cross  beam,  complete  the  output  data  base. 

Ten  strain  gage  conditioners  power  ten  bridges  and  continuously  sample  signal  out  from  them. 
These  conditioners  and  two  accelerometer  charge  amplifiers  are  located  at  the  center  of  the  crossed 
beams.  Line  power  to  these  conditioners  and  amplifiers,  and  output  data  from  them  transfer  through 
slip  rings. 

The  ten  data  channels  pass  into  ten  separate  filtering  circuits  that  rectify  and  eliminate  the  dc  bias. 
The  dc  bias  is  the  steady  component  of  total  signal  representing  the  centripetal  force  along  Fx  and  Mv. 
The  ac  component  of  the  original  signal  is  the  desired  sample.  An  ac  filter  attenuates  the  noise  signals 
above  4  Hz  The  rms-to-dc  chip  takes  the  absolute  value  of  this  signal,  multiplies  it  by  0.707,  then  con¬ 
verts  it  to  a  dc  result  representing  the  rms  value. 

Phase  Angle  Determination 

Two  accelerometers  provide  the  tank  motion  reference,  linear  and  angular  motions  are  sensed 
One  accelerometer  is  mounted  on  the  top  horizontal  ting,  outboard  side,  and  the  other,  on  the  bottom 
horizontal  ring,  inboard  side.  The  positive  axes  of  both  accelerometers  are  mounted  along  the  vertical 
in  opposite  directions  Subtraction  of  linear  components  result  in  a  reference  signal  proportional  to  the 
angular  acceleration,  a,,  about  the  y-axis.  The  phase  of  the  A/,  signal  relative  to  the  a,  signal  was 
then  determined  by  electronically  measuring  the  time  difference  between  the  peaks  of  the  two  signals 
In  practice,  the  A/,  gage  system  senses  the  moment  of  all  forces  and  couples  above  the  spherical  bear¬ 
ings  (Fig  5).  Thus  the  desired  signal  proportional  to  the  viscous  shear  couple  about  the  tank  center 
was  not  measured  directly,  which  somewhat  diminished  the  usefulness  of  the  phase-angle  measure¬ 
ments  as  a  method  of  estimating  energy  dissipation  rates. 

TEST  PLAN 

Tests  represented  conditions  with  water  or  glycerol  as  the  sloshing  liquid  in  tanks  located  at  32  in. 
radius  and  21  or  22  in  above  spin  center.  The  lower  support  ring  of  the  control  tank  was  weighted  with 
lead  equal  to  liquid  weight.  Beginning  with  the  desired  fill  fraction  (50,  75,  or  100%)  and  nutation 
angle  (1°,  2°,  or  3°),  a  constant  outer  table  spin  rate  (0,  15,  20,  25,  30,  60,  80,  92.5,  107  rpm)  was 
selected.  The  nutation  (inner)  table  was  rotated  from  10  to  100  rpm  and  at  2-rpm  intervals,  the  sen¬ 
sors  were  sampled.  Each  nutation  frequency  was  held  constant  for  a  time  sufficient  to  allow  transients 
to  decay.  Data  were  plotted  by  computer  after  completing  a  test.  This  readily  showed  a  liquid  slosh 
resonance  in  which  the  increased  force  and  moment  observed  in  the  liquid  lank  was  net  seen  in  the 
control.  Not  every  test  condition  produced  one  of  these  resonances.  Another  test  condition  was 
selected  and  the  test  repeated.  All  results  were  used  for  boundary  layer  shear  computations. 

Fifteen  of  the  81  tests  with  water  were  repeated  with  glycerol  solution  (with  viscosity  ten  times 
greater  than  water)  to  verify  selected  model  damping  parameters. 
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EXPERIMENTAL  RESULTS 

Figure  6  shows  a  typical  slosh-wave  resonance  excited  at  a  low  spin  rate  and  high  nutation  rate. 
Here,  F..  exhibits  the  classic  resonance  response.  We  note  that  its  moment  pair,  Mx ,  showed  similar 
shape.  The  Fx  and  M>  pair  were  mirror  images  of  these;  their  magnitudes  decreased  before  sharply  ris¬ 
ing.  This  was  affected  by  rotation  of  the  nodal  line  of  the  surface  wave  around  the  tank. 


LIQUID  - 

CONTROL  - 


lo  10  20  40  60  80  100 

X  NUTATION  RATE  (rpm) 

SPIN  ANGLE  =  20  RPM 
CONE  ANGLE  -  1  deg 
FILL  FRACTION  -  60% 

TANK  DIAMETER  -  12  inches 
RADIAL  LOCATION  -  32  Inches 
VERTICAL  LOCATION  -  21  inches 

Fig  6  -  Typical  low  spin-speed  result  showing 
slosh-wave  resonance 

Figure  7  shows  responses  excited  at  high  spin  rate  and  low  nutation  rate.  All  high  spin  rates  (ft 
^  60  rpm)  excited  a  structural  resonance  in  the  table.  Note  increasing  amplitude  as  a  function  of  9 
and  that  changing  viscosity  does  not  greatly  affect  amplitude. 

Calibration 

Voltage  outputs  of  the  strain  gages  were  calibrated  by  using  different  methods,  depending  upon 
table  resonance 

The  table  did  not  resonate  at  low  spin  rate  (fl  $  30  rpm).  In  this  spin  region,  the  dynamic  cali¬ 
bration  method  required  determining  the  expected  forces  and  moments  sensed  by  the  dynamometers 
according  to  the  principles  of  dynamics,  for  nutation  frequencies  well  below  the  slosh  resonances.  Cali¬ 
bration  curves  resulted  from  loading  the  empty  liquid  tank's  lower  horizontal  ring  with  lead  weight 
equivalent  to  50,  75,  and  100%  fill  fractions.  Then  expected  forces  and  moments  were  divided  by 
actual  voltage  outputs  at  specific  spin  and  nutation  rates  to  farm  calibration  constants. 

At  high  spin  rates  f O  >  ■')  rpm),  the  table’s  structural  resonance  effectively  increased  the  ampli¬ 
tude  of  the  table  motion  of  the  tank  location  in  an  unknown  way.  Therefore,  there  is  no  direct  way  to 
calibrate  the  signals.  Instead,  an  indirect  semianalytic  method  was  used  To  cancel  the  effect  of  the 
unknown  table  amplitudes,  the  voltage  representing  My  was  divided  by  the  voltage  representing  .  * 
(and  similarly  for  MJF  )  at  each  nutation  rate.  The  ratios  were  normalized  by  dividing  them  by  the 
corresponding  ratio  when  the  nutation  rate  equals  the  spin  rate,  X  -  Fl  at  this  condition,  the  tank  does 
not  oscillate  but  merely  spins  steadily  around  an  axis  inclined  at  9  to  the  vertical,  a  condition  for  which 
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SYMBOL  LIQUID  HALF-CONE  ANGLE 


•  WATER  1 


SPIN  RATE  =  60  RPM 
FILL  FRACTION  =  50% 

TANK  DIAMETER  =  12  INCHES 
RADIAL  LOCATION  =  32  INCHES 
(CONTROL  TANK  DATA  NOT  SHOWN) 


Fig  7  —  Typical  high  spin-speed  results 
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the  true  moments  and  forces  are  easily  calculated.  From  these  normalized  ratios,  along  with  the  known 
moments  and  forces  at  X  -  ft ;  the  moment  and  force  data  could  be  calibrated. 

Determining  Model  Parameters 

Slosh  resonances  were  excited  in  tests  driven  at  low  spin  rate  (ft  ^  30  rpm)  and  high  nutation 
rate  (X  >  80  rpm).  Figure  6  shows  a  typical  result  of  the  oscillating  force  magnitude.  The  nutation 
frequency  at  which  the  resonant  peak  occurs  is  used  to  determine  the  pendulum  length  /,.  (Only  the 
lower  of  the  two  slosh  resonances  mentioned  in  the  third  section  of  this  report  was  observed  because 
the  second  frequency  is  above  the  nutation  capability  of  the  table.)  The  height  of  the  resonant  peak 
force  above  the  equivalent  rigid-liquid  baseline  and  the  width  of  the  resonant  curve  are  sufficient  to 
determine  the  pendulum  mass  Ms  and  the  viscous  dashpot  C,  (3). 

The  magnitude  of  the  rigidly  attached  mass  A/0  equals  the  difference  between  the  total  liquid  mass 
and  the  pendulum  mass  (MT  -  Ms).  From  considerations  of  symmetry,  both  M0  and  the  pendulum 
hinge  point  must  be  located  at  the  center  of  the  tank,  independently  of  the  fill  level  (3],  Similar  data 
analyses  determine  the  parameters  of  the  pendulum  model  for  all  fill  levels,  spin  rates,  and  liquid 
viscosities.  The  results  for  A/0,  A/j,  and  /,  are  not  markedly  different  from  those  found  previously  for 
a  nonspinning  spherical  tank  (3)  with  the  exception  that  the  viscous  damping  is  three  to  five  times 
greater  for  the  spinning  tank.  A  nonlinear  liquid  response  at  resonance  complicated  data  reduction. 
The  nodal  line  of  the  surface  wave  rotated  around  the  tank.  Similar  behavior  has  been  observed  in 
nonspinning  spherical  tanks  (3J. 


Inertial  wave  resonances,  if  they  exist,  should  be  excited  in  tests  driven  at  large  spin  rates  (ft  ^ 
60  rpm)  and  small  nutation  rates  (x  <  40  rpm).  The  data  did  not,  however,  reveal  any  resonances  that 
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were  large  enough  to  note  in  the  moment  measurements.  Because  the  moments  were  measured  about 
an  axis  that  did  not  pass  through  the  center  of  the  tank,  the  data  had  to  be  analytically  adjusted  to  sub¬ 
tract  out  the  large  moment  of  the  oscillating  liquid  force  about  the  dynamometer  axis.  Even  though 
this  type  of  data  reduction  has  some  inherent  inaccuracy,  there  was  no  evidence  of  resonances.  After 
correcting  the  data  as  indicated,  the  moment  about  the  tank  center,  which  is  due  to  the  viscous  shear  at 
the  wall,  can  be  determined.  Figure  8  shows  a  typical  result.  The  straight  line  fitted  through  the  data  is 
the  type  of  correlation  expected  for  an  unsteady  viscous  boundary  layer,  in  which  the  bulk  of  the  liquid 
is  unperturbed  from  its  state  of  steady  spinning.  (The  data  shown  in  the  plot  do  not  exactly  fall  on  a 
straight  line,  even  allowing  for  the  rather  large  scatter;  but,  the  curvature  of  the  data  is  not  always  evi¬ 
dent  in  other  tests  and  in  some  tests  displayed  the  opposite  curvature.  A  straight  line  relation  implies 
that  the  viscous  shear  is  linearly  dependent  on  the  oscillatory  rotation  rate  of  the  tank  and  represents 
the  best  correlation  for  all  the  data.)  This  type  of  moment  data  is  sufficient  to  determine  the  magnitude 
of  the  viscous  dashpot  that  represents  the  viscous  boundary  layer  at  the  wall.  The  best  fit  correlation  is 

C  =  34  1  (v/ild2)12  in-lbf-s ,  (2) 

where,  0  <  <  1  is  the  fraction  of  the  wall  wetted  by  the  liquid,  FK  depends  on  the  fill  level.  The 

measured  viscous  shear  is  several  orders  of  magnitude  larger  than  would  be  predicted  by  a  laminar  flow 
analysis,  this  result  is  in  agreement  with  data  presented  by  Vanyo  [6,7].  The  energy  dissipation  rates 
are  also  correspondingly  larger. 
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Fig.  8  —  Wall  shear  moment  for  half-full  tank  at  60-rpm  spin  rate 


The  apparent  lack  of  an  inertial  wave  resonance  implies  that  viscous  coupling  of  the  liquid  bulk  to 
the  oscillatory  tank  rotation  (which  is  simulated  by  the  dashpot  C  in  the  model)  is  too  weak  to  cause  a 
noticeable  deviation  of  the  liquid  from  a  steady  state  spin.  This  observation  agrees  with  previous  quali¬ 
tative  tests  [8],  The  mechanical  model  also  predicts  than  an  inertial  wave  resonance  in  a  spherical  tank 
would  be  rather  small,  regardless  of  the  magnitude  of  the  dashpot  constant.  That  is,  the  dashpot  (liquid 
viscosity)  drives  the  liquid  motion  as  well  as  damps  it.  If  the  dashpot  constant  is  large,  the  driving 
moment  on  the  rotor  (that  is,  the  liquid)  would  be  large  but  so  would  the  damping.  Conversely,  if  the 
dashpot  constant  is  small,  the  damping  is  small  but  so  is  the  driving  moment.  Thus  the  resonance  will 
always  be  a  small  effect  in  comparison  to  the  large  resonances  driven  directly  by  rotation  motion  of 
nonspherical  tanks  [4], 

Not  all  inertial  wave  model  parameters  could  have  been  determined  by  test  data  even  if  an  inertial 
wave  resonance  had  been  observed.  The  resonant  peak  moment,  resonant  frequency,  and  resonant 
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curve  width  would  only  be  sufficient  to  determine  certain  ratios  of  model  parameters,  such  as 
lt/lx,  ljly,  kxCJDx ,  and  kyCy/Dy.  These  ratios,  however,  are  all  that  are  needed  to  compute  energy 
dissipation  rates  both  at  the  resonant  frequency  and  at  nutation  rates  far  removed  from  resonance. 
Thus  the  lack  of  definiteness  would  not  be  a  drawback  in  applications.  The  uncertainness  of  the  model 
parameters  does  point  out,  nonetheless,  the  need  for  a  complete  theory  of  liquid  motions  in  a  spinning, 
nutating  spherical  tank. 

SCALING  RELATIONS 

One  of  the  advantages  of  an  equivalent  mechanical  model  in  interpreting  and  generalizing  labora¬ 
tory  data  is  that  the  model  itself  can  be  used  to  predict  on-orbit  conditions;  otherwise,  the  test  data 
would  have  to  be  scaled  directly.  The  way  the  model  parameters  depend  on  spin  rate,  fill  level,  viscos¬ 
ity,  and  other  tank  and  liquid  properties  can  be  established  more  reliably  than  can  the  scaling  of  gross 
measurements  of  energy  dissipation. 

By  considering  the  physical  variables  that  describe  liquid  motion  in  a  spinning,  nutating  spherical 
tank,  a  set  of  nondimensional  scaling  parameters  for  the  model  has  been  formulated  by  dimensional 
analysis  [9],  Test  data  were  used  to  determine  the  importance  of  the  various  non-dimensional  parame¬ 
ters  and  to  quantify  the  empirical  constants.  The  scaling  relations  for  the  pendulum,  or  slosh  wave, 
model  are: 

MT , 

1,/d ,  and 

CJM.l'Q.  (vV^fc//)01*. 

where  ge//  -  Lg2  +  (Jr0ft2)2ll/2.  Numerical  values  of  these  ratios  determined  from  the  tests  can  be 
used  in  the  model  to  compute  full-scale,  on-orbit  energy  dissipation  rates.  For  the  rotor,  or  inertial 
wave  model,  only  the  boundary  layer  dashpot  constant  was  found  to  be  important.  This  is  scaled  as  fol¬ 
lows: 

c/prf5flf»WftrfJ)l/I. 

The  other  parameters  of  the  rotor  models  should  scale  as  shown: 

/.//,, 

/A. 

kX,  ky  , 

ex,  e,,  and 
D/pdSfUv/nd2)* . 

However,  since  no  inertial  wave  resonances  were  found,  numerical  values  of  these  parameters  could 
not  be  established. 

SUMMARY 

This  report  describes  a  laboratory  based  method  to  estimate  energy  dissipation  rates  for  prolate 
spinning  spacecraft.  A  forced  motion  drive  table  was  designed  and  constructed  to  spin  and  nutate  scale 
model  propellant  tanks.  Because  the  nutation  frequency  can  be  adjusted  independently  of  the  spin  rate, 
the  table  is  able  to  excite  all  important  liquid  resonance  frequencies  and  measure  forces  and  moments 
throughout  the  range  of  interest. 
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The  experimental  data  were  used  to  determine  the  parameters  of  an  empirical  mechanical  model 
of  the  liquid  motion  dynamics.  In  the  model,  a  pendulum  represents  the  free  surface  modes  of  oscilla¬ 
tion  (sloshing)  and  a  rotor  represents  the  inertial  wave  modes  of  oscillation.  The  rotor  is  connected  to 
the  tank  by  a  dashpot  that  simulates  the  forcing  of  inertial  waves  by  the  boundary  layer  shear  at  the 
tank  wall.  Other  dashpots  are  used  to  simulate  the  viscous  damping  of  the  slosh  and  inertial  wave 
modes. 

The  test  results  and  model  determinations  have  shown  that  the  sloshing  modes  in  a  spinning 
spherical  tank  are  not  significantly  different  from  those  determined  previously  for  nonspinning  tanks, 
with  the  exception  that  the  damping  is  increased  by  a  factor  of  three  to  five.  The  tests  also  show  that 
no  inertial  wave  mode  is  excited  to  a  significant  amplitude  by  the  boundary  layer  shear  in  a  spinning, 
nutating  spherical  tank.  However,  the  boundary  layer  effect  and  the  corresponding  energy  dissipation 
are  much  larger  than  would  be  anticipated  on  the  basis  of  a  laminar  flow  analysis. 
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Appendix 

ANALYSIS  OF  MODELS 


COORDINATE  SYSTEMS 

Figure  A  1(a)  shows  that  the  equilibrium  (nonnutating)  position  of  the  tank  center  is 
X'  =  A'q,  Y'  =■  0,  Z'=  Zq,  referred  to  a  rotating  X\  Y\  Z'  system  located  at  the  pivot  point  of  the 
forced  motion  drive  system.  The  equations  of  the  model  are  developed  in  a  body  fixed  set  of  axes 
x,  >\  z,  which  for  clarity  are  referred  to  an  X ,  Y,  Z  rotating  system  similar  to  X\  Y\  Z'  but  located  at 
the  tank  center.  The  angular  rates  about  the  x,  y,  z  axes  to  the  first  order  in  the  half-cone  angle  6  are 

wx  *  -  (H  -  k)6e,KI  (Ala) 

(x)y  =  ~ii.il  -  \)6e ,  and  (Alb) 

id2  =  fl  ,  (Ale) 

where  only  the  real  parts  have  physical  meaning. 

PENDULUM  MODEL 

Figure  Al(b)  shows  that  the  linearized  equations  of  motion  of  the  pendulum  are  expressed  in 
terms  of  small  angular  displacements  a  and  d>.  The  displacements  are  referred  to  the  equilibrium  posi¬ 
tion  yo  of  the  pendulum  as  determined  by  the  combined  action  of  gravity  g  and  centripetal  acceleration: 

tan  >0  =  g/ln^Afo  + /,  cos  y0))  • 

At  equilibrium,  the  pendulum  lies  in  the  x,  z  plane. 
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By  considering  moments  about  the  hinge  point,  the  equations  of  motion  for  a  and  <p  can  be 
shown  to  be 

4  la  +  n a  a  +  2fl  iji  sin  yo  cos  y0l 

-  -  0o  sin  y0  l#  +  (Z o  -  4  sin  yo)  (fi  ~  X)2)ea' ,  and  (A3a) 

4  cos  y0l«//  +  nl  ill  -  2fld  tan  y0) 

"  +  (Z0  -  4  sin  Vo)  (H  -  X)2]e,A' ,  (A3b) 

where 

fj.2  -  n:  (cos2y0  +  (AV 4  cos  y0)l ,  and  (A4a) 

nZ-Cl2(X  o/4cosy0).  (A4b) 

The  effect  of  damping  has  been  neglected  in  this  derivation,  but  as  customary  it  can  be  added  later  in 
the  resonance  terms.  Letting  ao  and  i//0  be  the  complex  steady-state  amplitudes  of  a  and  ip,  the  solu¬ 
tions  of  Eqs.  (A3)  including  damping  are 

ac  -  0[(jr2  ~  X2)i?0  —  2 fl\A0  sin  yol/Go»  and  (A5a) 

ipo  -  0(2fixBo  (an  yo  -  (mJ  ~  X2)Ao/c os  y0)/6o<  (A5b) 


where 


where 


A0  -  (g/ls)  +  [(Zq/Is)  -  sin  y0](ft  -  X)2, 

B0  -  ((Zo/4)  +  cos  yol(X2  -  H2)  -  A0  sin  y0,  and 
G0=  (X,2-  X2  +  2/XX,£,,)(X22  -  X2  +  2 /XX2£j2). 


(A6b) 

(A6c) 


The  damping  ratios  are 


£ji  "  \ ,  and  (A7a) 

£32  ”  Cj/2A/j42X2,  (A7b) 

and  the  pendulum  natural  frequencies  are 

(x  12/ f! )2  “  (A’(/4  cos  y0)  +  y(l  +3  sin2y0) 

T  l4(Ar(/4  cos  y0)  sin2  y0  +  -*-(1  +  3  sin2  y0)2],/2.  (A8) 

The  X|  and  X2  cannot  be  written  in  terms  of  an  "effective"  gravity,  say  [#02  +  (A'  fl2)2]1/2,  although  the 
error  involved  in  doing  so  is  not  large.  In  orbit  where  g  =  0,  Eq.  (A8)  reduces  to 

X,  -  n(AV4),/2  and  (A9a) 

X2=  fill  +  (AV4)],/2.  (A9b) 


(A9a) 


X2=  fill  +  (AV4)],/2.  (A9b) 

By  comparing  these  expressions  to  Eqs.  (A4),  X|  is  associated  with  a  resonances  (transverse  motion) 
and  X 2  is  associated  with  ip  resonances  (circumferential  motion). 

The  components  of  the  liquid  force  exerted  on  the  tank  can  be  written  in  terms  of  a  and  1//  For 
example. 


Fx  =  -Mr  [A0  +  cos  y0)  O2  -  MTg8e'" 

-A/t-IZo  -  ( MJ Mt)Is  sin  y0I  (f l  -  a)2  ve:K! 

—  ( Msls[(d  —  afl 2)  sin  yo  +  2fli//  cos  yo))  • 


(A  10) 


»>  *  s  »• 


r-v 

V  *  -  v  •*. 

"I 

*■ 


-.-A 

-X.-A-s 
a/  A*! 
«  ■ «. « 


LJI 

vrji 
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The  term  in  braces  represents  the  effect  of  pendulum  motion  (sloshing).  The  other  terms  are  merely 
the  rigid  body  response  of  the  total  liquid  mass,  MT  -  M0  +  Ms,  concentrated  at  the  center  of  mass  of 
the  liquid. 

The  energy  dissipated  per  nutation  cycle  of  the  liquid  motion  is 

-  nkCs(  |c*ol2  +  l»//olJ)  ■  (All) 

ROTOR  MODEL 

The  angular  momentum  of  the  steadily  rotating  liquid  is  I:Cl,  where  /.  is  the  moment  of  inertia 
of  the  liquid  about  the  z  axis.  When  inertial  oscillations  are  excited,  additional  angular  momentum  is 
created.  This  additional  momentum  about  the  tank  center  can  be  expressed  as  lx  hx  and  Iy  hy,  where 
hx  and  /),  are  momentum  averaged  angular  rates  of  liquid.  (In  the  linear  approximation,  the  additional 
momentum  about  the  z  axis  is  neglected.)  lx  and  /v  may  be  thought  of  as  moments  of  inertia,  but 
they  do  not  have  any  rigid-body-like  t elation  to  /.  and  the  geometrical  shape  of  the  tank.  If  lx  and  Iy 
can  be  approximated  as  time  independent,  the  fluid  motion  can  be  modeled  as  a  rotor  at  the  center  of 
the  tank  The  equations  of  motion  for  the  rotor  are: 

Ix  rix  +  Dx  hx  +  (/.  -  Iy)  Cl  hy  -  Tx  and  (A12a) 

Iy  Hy  +  Dy  W,  ~  (/.  -  /,  )  0  Hy  “  Ty  ,  (A^) 

where  Tx.  7",  are  the  moments  applied  to  the  liquid  by  the  viscous  boundary  layer  at  the  tank  wall,  and 
Dx,  Dy  are  the  x\  y  components  of  the  dashpot  shown  in  Fig.  1. 

The  viscous  moments  are  modeled  as 

Tx  —  C^O  +  /)  [<ux  -  kxhx  —  iex  ( Cl/\ )  ( qjv  —  ny ) )  and  (A13a) 

Ty  »  Cv(l  +  /)  [w,  -  kyhy  +  iey  ( Cl/k )  (wx  -  fcx/r*)],  (A13b) 

where  Cx ,  Cy  are  the  components  of  the  dashpot  C.  The  term  (1  +  /)  represents  the  fact  that  an  oscil¬ 
lating  boundary  layer  has  components  of  shear  that  are  in  phase  with  both  the  wall  velocity  and  the  wall 
acceleration*.  Since  the  liquid  "free-stream"  velocity  at  the  wall  may  be  greater  or  less  than  that  indi¬ 
cated  by  the  average  rotation  rates  hx  and  hy,  empirical  factors  kx  and  ky  are  introduced.  Considering 
the  smallness  of  the  viscous  coupling,  an  inertial  mode  that  involves  a  lot  of  liquid  motion  in  the  inte¬ 
rior  would  be  one  that  needs  only  a  small  velocity  at  the  wall,  that  is,  the  magnitude  of  kx  and  ky 
would  be  much  less  than  one.  The  empirical  parameters  ex  and  ey  bring  in  the  Ekman  like  skewing  of 
a  rotating  boundary  layer.  In  the  following  development,  ex  and  ey  are  set  equal  to  zero  for  conveni¬ 
ence,  but  this  is  not  a  necessity. 

The  steady  state  amplitudes  of  the  oscillating  rotation  can  be  found  from  Eqs.  (A12)  and  then  put 
into  Eqs.  (A13)  to  predict  the  moments  exerted  on  the  tank.  Assuming  that  I kx  Cx  I  <<  1  and 
I ky  Cy  I  <<  1  (although  neither  Cx  nor  Cy  need  be  small),  the  results  are: 

Tx  —  -Cx(i  +  /)  (fi  —  k)  9eiK>  x  ( A 1 4) 


and  a  similar  expression  for  Ty  with  the  x  and  y  subscripts  interchanged  and  6  replaced  by  i6.  The 
symbol  Ho  represents  the  oscillation  effect: 

*H.  Schlichting,  Boundary  Layer  Theory  (McGraw-Hill,  New  York,  1979),  7ih  ed. 


ZEDD  AND  DODGE 

l,  i  2  ,  kxCx  *,C, 

//°  -  1  -  —  -2 /  C,  +  {  v  +  -ry- —  +  rj  ,  where 

*  *■  ‘x  ^  n 

”  DJ2lx\„ ,  and 
C,  -  A/2/vxn 

are  the  damping  ratios  of  the  inertial  oscillation  mode.  The  natural  frequency  of  the  mode  is: 


n  -  A_^l|  U — L 


(A15) 


(A16a) 
( A 1 6  b ) 


( A 1 7) 


( LHX )  >  1  and  (/,//>.)  >  1  for  stability . 


When  the  nutation  frequency  is  not  close  to  the  natural  frequency,  the  effect  of  the  rotor  is  small 
and  the  moment  on  the  wall  reduces,  for  example,  to 


Tx  -*  -Cv(l  +  i)  ( n  -  A.)  8elK' 


(A18) 


Equation  (A18)  simulates  the  viscous  boundary  layer  shear.  At  resonance,  where  A  =  A„,  the  wall 
moment  works  out  to  be 


7,  -  —  C,  ( 1  +  /)  (a  -  A)  0e‘K'  Jl  + 


(1+1)  1 1  c>-  I  n  I  f  /•  -  /> 


U  (A  1 8a) 


where 


K  _  _ _ ( k  x  CJ  /A  A  „ ) _ _ 

0”  (Cx  +  {>  +  (kxC,//xK„)  +  (kyCy/ly\„)  '  (A18b) 

Since  {x,  £v,  Cx,  and  C,  all  represent  viscous  effects,  it  can  be  seen  from  K0  that  the  resonance  is 
bounded  (and  small)  whether  £x,  Cx,  and  Cv  all  approach  zero  or  all  approach  infinity.  Wall  shear 
is  unaffected  by  the  inertial  oscillation  if  Cx  -  C,  and  fx  -=  ly  since  the  term  multiplying  A0  is  zero 

The  energy  dissipated  by  the  inertial  oscillations  is 

E,  -  rr  (0  -  A)  9  (|7*J  +  |  Ty,  I  )/A  ,  (A19) 

where  Txr  is  the  real  part  of  the  amplitude  of  Tx  and  Ty,  is  the  imaginary  part  of  the  amplitude  of  Ty. 


